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A  Jj  be  an  nxn  matrJ:'’  of  re-ii  ji'rnbHr.i.,  The  travel  1  mg  ,".aies~ 

II  ^Jii 

ruaii  problem  aaks  for  an  ac^rtie  ixirmutation  <1^  i  '  of  the  mbe;  ■  . 
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a  rainimum  .  If  ■  a.^  for  atl  J  and  .5,  the  problem  is  saici  to  be  sjm- 
iagtrli-;o  otherwise,  it  Is  nono^mmetrlc  In  t  a&a  the  distancea  ^  arfi  -om- 
puted  between  points  in  the  iluclidaon  p  amy  the  protlem  ts  du-sl  Idean..  n  r.'n 
clidean  prob  em  is,  of  course,  aymmetri*' . 

In  the  nonsvmmetrl-;  -iaBe  the  aotu  ion  involves  finding  one  ou)  oi’’ 

(ii  D  '  possible  perroutationat,  and  in  the  u.nmae -  ri  -  ase  it  j;.nwaive&  finding 
:>rxe  euf  o.t  r  ( n  I):  possible  pe-rmut-atxons  H&ni  e  .  >;ompi.att"  enuDieraTion  and 
evaluation  of  all  the  possible  permutationi  proviueM  a  theorettoally  satia- 
t'a>--tory  solu’ ion  to  the  problemi,.  What  as  actiuilly  desired-  however,  la  a  com- 
putalionally  pra  hi.ai  method  of  finding  the  optiuuim..  di  en  methods  that  find 
•'gond*'  permutations  would  te  of  practical  interest  as  appi-oximata  so'lutionso 
The  nam*s  "teaveli xng  oalesraar."  is  appll cable  when  one  interprats  the  S'um 
!'l)  as  tht  total  distar.f.t  that  must  be  travelled  by  a  salesman  who  must  visit. 


;ach  of  n  ‘-itiss  exaoily  once  before  rntuning  home-  The  acyclic  n-  xirement 
prevents  a  solution  involving  several  disconmnted  loops.  The  travelling 
salesman  intf'rpretation  usually  results  In  a  symmatrir*  problem.  However,  if 
we  think  of  a  school  bus  that  has  a  given  number  of  corners  at  which  to  stop 
to  pi'.ik  up  c.bildren  in  a  city  which  lias  a  substantial  number  of  one  way  streets, 
3  nonsymmetric  problem  results,,  since  the  distam  e  from  a  t.o  b  may  well  be 
iJlfferant  from  the  distance  from  b  to  a  Still  anolhar  nonsyaffiatri*’'  iiiterpre 
bation  Is  that  of  a  machine  tool  that  parforma  a  gi  /a  set  of  jobs  repaTi lively- 


The  diatance  from  job  a  to  job  b  la  the  setup  »308t,  which  «san  be  different 
from  the  setup  cost  going  f-'om  Job  b  to  Job  a..  For  instance,  tn  a  paint 
mixer.  It  Is  easy  to  go  from  a  iighit  to  a  darker  soJoFm  but  diffl^suit  to  go 
in  the  reverat  direfjtlon  There  are  numerous  other  practissal  Interprefca- 
iiona  of  this  seemingly  frivolous  probJeiUv 

In  thie  papier  we  shaiJ.  discuss  a  method ,  suitable  for  ele4tronic  <?om- 
puters*  that  hJis  proved  capable  of  quickly  obtaining  solutions  for  problems 
having  about  oO  cities  or  less  in  symmetric  and  some  nonsymmetric  problems i. 

In  principle  the  method  can  be  used  for  any  size  problem.  Although  the  code 
does  not  guarantee  finding  the  optimum  tour.  It  can  be  used  over  and  over 
several  times  and  in  various  ways  to  get  a  probabilistic  idea  of  how  good 
the  beat  answer  found  is  ri^latSve  to  the  set  of  observed  answers.  Also  the 
chs'-'king  procedure  of  Danfc_sig,  Fulkerson,  and  Johnson  l2]  can  be  appJlad  to 
check  on  the  optimality  of  the  observed  result. 

V/e  call  the  method  a  heuristic  one  because  ^a)  the  code  for  it  contains 
probabilistlr-  elements  so  that  its  performance  varies  each  time  it  is  run; 
fb)  in  certain  cases  It  can  be  proved  that  it  has  positive  probability  of 
producing  the  optimum  answer,  and  our  experience  leads  us  to  believe  that 
it  will  always  do  so,  and  ^v)  as  the  result  of  initial  calculations  partial 
answers  and  sub-problems  are  obtainewl  sc-  that  later  calculations  depend  upon 
the  results  of  early  'alcuiations,  i.e.. ,  the  process  is  one  of  "learning  from 
experience,"  Nevertheless,  the  code  is  an  algorithm  in  the  sense  that  it 
terminates  after  a  finite  number  of  steps  and  hae  been  run  on  an  electronic 
computer.  This  application  is  an  example  of  artificial  intelligence,  that  is, 
the  use  of  a  computer  to  solve  problems,  the  solution  of  which  by  human  beings 
would  be  regarded  as  intalligant  acts.  Humans  are  not  good  at  solving  trav¬ 
elling  salesman  problems  because  of  limited  arithmetic  abilities.  Hence  our 


codt  oo'iii  not  behaviur  of  h'SKianit  and  s-j  not.  iieuristie  m 


..CASp. 

A  gooA  ntjjianary  of  nhe  hi-atory  v'f  the  proolem  np  to  ;Jbi»  ie  pre'sented  'o 
Floo  ■  [*jj  ''ore  re^'entty.  fucke'*  [o']  afj-i  kaii..74j>5  [kj  have  given  Jjtteger  pro 
gi aiiiis' r'g  icirmo'ift  ions  of  t.ht  pr«>i>lL'io  and  sorao  ^offiputationaJ  exp<«rien.r«^  has* 
bmP-  gained  ul^h  inera.  The  largest  problm  -jioived  so  far  ir*  the  literature 
Is  t.ijtr  42  ci'y  pi’obioja  of  Lantzigs  t'jjfc^rson.!  and  Johnson  l2J..  unr  algorStto, 
also  solvyd  the  sama  42  city  probleta  in  4  i  minutes  on  a  Bendix  G-2U  »wza- 
pu-.r,.n,  4  mt‘h  more  difl'l-.ult  probleu)  Involving  H'7  cities  was  solved  .in  3% 

fn.inu'caa  Utter  experlonco  wj.'h  these  and  sraitlii^r  problacos  in  reported  on 
later  in  ‘>He  papej', 

i  1  Section  2  w»  aeacrib^-.  the  eimpie  idea  needed  for  -.he  basic  at@p  in 

t.hc:  progr.iEi  Iha  initial  code  wiiach  does  not  nave  btu.?."  in  iaarnliLg  is  dij>- 

oussad  in  tk'ction  The  final  *  <xle  togethar  witli  its  loai*ning  as;iect  is 

JlKOosad  Jn  Section  h  Jd  hcsttlotj  ‘j  the  pt oU_bi.l5;sil'‘  '.vihob^  we 

*jr.o  ■  '.■p  Ui.<  ansi/er**  invid,  aii-i  Li  -.H^dtion  L  we  prove  that  the  opt.baal 

tour  has  rx:>as.MvK  probabi J  ity  of  o-jtng  :ho-ien  with  the  a"gor.\t,iKjit 

"1]S  3.\;:JC  STik  UF  Tits  A-UloBJT«M 

In  the  Auclldoan  travelling  salesman  p'*cibJeni  H.  can  be  proved  that  an 
optiaal  tour  wil.  aevor-  c-roi;-.  itsi.lf  'Hd  follows  from  the  liuclidean 
theorem  that  the  s'jm  of  two  n.cdee  of  a  i^langlui  I--;  greater  than  the  third 

side  However,  wbt^n  ci<le«  on  a  map  art-  u«ea  the  curvature  of  the  earth, 

the  -^xlofenLn  of  juuiuitaiJis  having  paase-n  and  tunae.!K,  and  the  existence  of 
Likert/  oci^anc..  aixl  other  natural  'oarrlera,  negate  the  above  liucliaean  result. 

Of  course,  nothing  like  It  need  be  true  for  no.nsyimnetri ?  problems.  In  any 
casfu  if  ar'  optima  I  tour  can  be  found  by  any  means,  it  will  have  the  desired 
propvrtias  without  spatifically  stating  reatr let, ions  needed  to  insure  getting 
t.h'.m 


4 


For  thi.8  reason,  ihe  only  specSi’ir  reetrii';t4.an  on  t.h&  pBrmutations  w«» 
-:onatrut:.t  is  that  they  be  aeyclif;.  Wt  have  devised  an  inductive  method  for 
•’ont.txU'cting  a-ycilc  permutations ..  Vh(4  method  ts  desrribed  in  the  followini.' 


series  of  i'iteps: 

t  Chooae  any  two  cities  and  list  them  arbitrarily  to  form  an  aeyclia 
pemuitation  i.,)  of  length  two  . 

2.  Assume  that  a  permutation  fi.,  i^  " '  ■  ijj)  oitiesK  where 

2  <  k  <  n  has  been  constructed-  Choose  one  of  the  remaining  cities^  call  it. 
city  h,  For  J  riu:.ning  from  1  to  n  compute  the  quantities 


h4.  , 


r‘i*i 


whei'ff  Vi'  define  i_  .  be  i.,  when  Jj  -  n,. 

4 1 

3,.  Lai  ba  any  value  of  j  su  h  that  .Is  a  minimum  of  the  q»iantit4.eJ3 

computed  in  2 

# 

4«  he  label  as  for  J  '  j  «  v .  »  n  and  label  h  as  » 

5  rife,  t-hua  have  'onstructed  a  permutation  (i,  ,  c.c  1  )  of  k»l  cities. 

If  k»l  “  n  stop;  otherwise  replace  k  by  k  J.  and  retura  to  step  2^ 

Briefly,  the  method  <>ons3?t3  of  starting  with  any  pair  of  cities  as  a 
permutation  of  ’ength  2.  then  inserting  a  third  city  in  such  a  way  as  to  mini¬ 
mize  the  length  of  the  result.lng  tour  on  three  cities;  then  inserting  a  fourth 
city  in  su-’h  a  way  as  to  minimize  the  resulting  tour  on  four  cities,  etc. 

This  is  our  heuristir  rule  for  constructing  acyclic  permutations. 

The  to'ir  resulting  ff*om  this  afsyclic  permutation  may  or  may  not  be  the 
optimum  one.  In  fact,  there  are  a  whole  set  of  possible  tours  of  various 
lengths  that  can  be  so  generated,  depending  upon  the  order  In  which  new  cities 
are  introduced  In  Section  6  we  present  a  proof  of  the  fact  that,  In  the 
Sue  t  dean  "asea  there  always  exists  an  order  in  which  to  inttoduce  the  citie.« 


Jto  'S'hi?.'  thi  aDov&  heuriHti  i-aXi--  wj.!'  produ  «  thtf  opfcima.1  lo'ir^  W«  have  nearly 
cii'  .;  .5i<ted  \ri  ijrif  what  vt  Ivlieve  to  be  the  optinial  tour  even  Sri 

the  nor,  &a  ikiii-Lin  -a'l?;!  and  jjo  we  ^•onje’-tiicre  that  the  same  result  hoIWa  tor 
ihtm-  HoWf»v6ir  ve  do  not  tiave  a  proof  of  t.hfe  fact  at  the  present  time, 

'ife.  f  iv?  cJty  probWs  of  Figure  i  wilt  tielp  ni’oatrate  the  method  as  we.tl 
Ld  ]ts  dii'fi  uities.  Here  the  optJmai  tour  is  12345.  which  has  Jengtb  148. 

The  orJy  other  ^-our  -cionstruotad  by  our  algorithm  Is  12453 -•  whidi  has  lsngt.h  152. 


F‘<gwrt' 

'a'o  rar.  25  rrlr-lf.  of  the  atgorUruD,  more  fully  described  in  Section  3»  using 
random  orders  for  introducing  new  Cities,  and  found  that  15  of  t-hem  resulted 
In  the  opijiial  ♦‘oiu  .  while  tb^  other  lO  produced  the  suboptiinal  tour.  Thus, 
tor  this  prob'em,  there  tfl  empirtcaj;  proDabil?t.y  of  6  of  finding  the  opti¬ 
mum  tour  by  random  procedure  . 

n-K-,  MPST  AJ-.GC-fiITHM.  CUtE  ' 

Tht  result  of  the  algorithm  of  the  previous  section  depends  on  the  order 
111  wni  h  lilt'.  art  intcodu«ied  in  step  2.  So  that  final  results  ware  nox 
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biastci  v-e  first  ;;rrang9d  the  -lities  in  a  random  ordffir  list  and  seje^ted 
'  it.lt  £J  from  5  his  Usu.  when  needed  for  step  2.,  W«  did  this  a  n»jmb@r  of  t")rif'.3. 
ijrlrtl int;  'ach  Jour  and  it, a  distance  earb  ■^.imso  We  (’,aXi‘  the  first  algor!'''-., 
t.ode  .  Por  anatanee;;  in  the  lo-  cJty  probiem  atudird  by  L.  L.  Baraiihei  [11, 
who8t'  data  is  given  in  ;ixhibit  J  and  whose  optimal  tour  is  g.«ven  in  jidchibit  3,. 
the  distrsbution  of  ••jorapletion  times  is  given  ,iji  Figure  2-  Note  that  the 
Wf!pir3v''a'>  probability  of  getting  the  optimum  tour  fwhose  length  is  3'78  miles) 
i*=  16.  Note  also  that  the  distribution  is  multi-modal  and  has  gap,s.. 

Sonifr  other  experimental  results  are  displayed  in  Figure  y., 


1  No.  of  Cities 

J 

1 

1 

j  iSmplricai  Prob.  of 

'  getting  optimum 

j  with  Code  i . 

] 

o6j  16  uOl  ,ou45 

Figure  3 


The  probability  rstimates  for  the  5..  10*  33*  and  5?  were  based  on  runs  of  100 
tv,a*-ih  and  that  of  the  a?  vxty  problem  on  a  run  of  225-  7’he  shortest  schedule 
observed  in  tne  5  7  cj.ty  case  was  331  miles  longer  than  what  we  believe  to  be 
the  opllmiBn.  Wt*  did  not  feel  that  the  erapirtnal  probability  of  getting  the 
host  answer  la  the  57  oity  ^sase  was  high  enough  to  continue  somputations  with 
Code  1  until  one  was  observed-  The  time  to  construct  one  tour  in  the  57  city 
•■•asf-  wax*  about  ,j0  seconds,  and  the  time  for  the  other  problems  proportionately 

From  the  experience  we  obtained.,  particularly  on  the  33  and  42  sity  prob- 
jg-ms,  w.'*  found  that  the  sbort.cpt  tours  so  produce!  tended  to  agree  pretty  wall 
aioind  i.h©  periphery  of  the  tour  whom  it  tended  to  be  convex*  but  did  nol 
at  all  well  in  the  ''canter"  of  the  problem  where  the  optimal  tour  waa 
necessarily  quite  non  convex.  Hence,  we  built  into  Code  1  the  added 


f  ii  '  •  i  j,  ot'  iig  tt.i  oi5^'ci-f„v  v^ul  .iAo  b.'vvi-  ‘i  liH'  ;.rrv .'i'afii  /i-j/i 

<m  t:; -ill  fi'.- »•• 

■.'■.'i'  Invj  ai'i  'fi  in  uir  haV'--  1  .■  a  ’n  wb}<  f;  the. 


l.r.e  i  .-'via  a  <-o).vex  par^-  of  ite  toui'  'rlu^t  waa  Irtiverst'd  by 

wo^v.  '•1  ■''V  V-;  01  f  ^oura  prc'ducsd  ‘n  tit  Snirial  eiinh  n.f  rbe  pri  blea-  Pcf 

l')<y  wo  riaci-  i  ''-ul'”  I'iopi  i.  '•€  :  aijd  >o'iv°d  f.(v>  tv*'}  i?uF-p‘r-oblt*i!e 

aa  pa.*H '■  *.  i  •:■  if'  '  hif  xny  W6  "fa  ■  tor  Ml"  t'H'  t)  .I'.v  probleiii  into  ail  ^ 

an  n  C  pr'.i  *.n  lnjritie  .anl  o-f.ur  jpi  ha  h  i.tfc.>hj.ti3.s  ,  Ihe 

"-■  •  ■.p"  pr  -  ‘-a,  fci'-ic  a  ■onv^'x  ''or  ncar-iy  '  prob.aii  ..i-  t?s»y 

i-:*  '•  :/i  'i.  iv.  i  aan  usuaJl;,  '"Uiultf  :  ;n  probabi.iJ.ly  of  1  or  greatsiir  io 

+  h-..  ■■.T';'')-  -.p,'.  -y  .  {■  rb'iHnsitir-n  ■•■•f  r'-.e  shcr-^eat.  Tj.n:  .  Thuis.*  at  wii,; 

0*  "Xr'  t..  'I*.;  in  ion  ”■  Vfv  wefs  virtually  ■.:3’-'air,  o,f  obialniny,  th*;  rarrev 

.ij-i'iW'.'-  o  .  Mh  -  or -bj  fME  T!;e  ri~6  profa-r'iii  ran  be  lonsicidred  ae  a 

pf '->b '  on  j’  ouai  c  to  fuftfer  f.ivt.or  » no  af  the-  saaie  kind  TFa  only 

ni  sf  e  :;c',  utioriii  to  ifto  u»rjd8  >vf  dibprcblftms  so  •lefin®  i  ir*  that 

thr  &:?:<;•'  i.n  Figu-f'e'  A,  the  T-oirmori  3 ink,  la  ahc-ioi  daubed  from  i,  t,o  1,, ! 

h.  c 

i:i\x»>  e  i  > .  .i-.i-ruMi  on.ce  ;n  the  uclufeiona  to  both  aubprcbiecif!,  Ii,  htw  ^ai-e  ^ 
th^-  ■ic-r;-y..w.-.(,r3.'  prc»bif«ii,  there  la  the  fiir'ther  ra^ibirarai^nt  Inat  ccmmon 


UnK,  b!'  (rweryea  cnoi>  in  owf  ■dl^e*rt4on  in  the,'  solution  to  the  firat  aub- 
i>''i  b.l  .-ca.,  .v.d  on;>'  in  tho  opj-xi^ite?  .iirti'=.tion  in  the  aotution  i<?  the  roi&ple- 
nc’ntary  i-ubpi-vb  'r-a. 

uur  t  xi'erat-'iiji..  as  exhibited  in  figure  i,  with  t-he  42  and  33  Jlty  prob- 
l?£riy  wap  that  wt  could  obtain  the  optimal  solution  to  ea^sh  problem  using  aXI 
t.faf  oitie.a  Bat  the  empirical  probability  of  finding  thera  in  this  way  was 
quit-?  sraalK  We  then  verified  probabllistical.^y  (see  Sesstion  5)  that  the  baat 
tour  was  optima*.,  by  factoring  them  into  sub  tours  and  showing  that  ae  sub-- 
problems;,  theue  subtoui's  !iad  empirical  probability  of  1  or  greater  of  being 
pi..ked  up  by  Cole  1 

Ooi*  rifeoct  Idea  was  to  have  the  vomputer  do  the  factoring  of  the  probl.em, 
arc  ■/?:  dfOfribf  how  it  did  that  In  the  next  aertion, 

4  c  t;-e:  ]  me  alg-jHIThm  .  joole  £ 

JUT  au'  cess  in  factoring  problems  aft-er  making  some  initial  runs  with 
Godt-  an  bough  aatlsfaatory,  was  t)me  cionsiiming.  and  led  us  to  incorporate 
tht  fafitcr.'jig  p’-oc&dure  into  the  code..  The  result  is  Code  2  to  be  described 
nei.t  !■<  g  a  -odt.  that  has  ’earning  aspects  an  that  the  results  of  earJy 
coi'putat:  ons  cf  the  program  define  the  subprobleos  chosen  to  be  worked  later 
in  m^-rp  detaij*.  Of  rourse.  both  eoi^rert  and  incoiToct  learning  are  possible 
Tht  a  djo?  1';  bni  wp  thus  developed  ie  one  that  stops  after  a  finite  number  of 
8tt  p6  a.u-  prints  out  a  tour  that  may  or  may  not  be  optimal.  Our  experience 
abnwj  .hfW  it  ha*'  positive  fflrobabillty  of  finding  the  optimal  tour  which  de- 
r?rid*:  ’itx'ri  the  ;<tze  of  the  problem  being  worked  and  on  the  number  of  random 
tr-.-3S  th;  iiiachJii':'  3fi  permitted  to  make  ‘'ioS.,  the  experience  it  is  allowed  to 
ha’  e)  before  bebig  requixed  to  define  subproblems  to  be  worked  on,  To  describe 
jutit  one  numerical  result  with  Code  2  we  found  that  in  the  42  city  problem  the 


algoritljm  found  the  correet  answer  after  making  5  cuts,  l^soi  it  defined 
6  nubproblems  ^  The  j'uts  -ran  be  observed  irx  iixhibit  D„  Other  results  wi:n  be 
discussed  iater. 

Vie  di?3f.ribe  in  inoni  dete-lJ  the  a,'goi  itlm  for  Code  2  whish  includes  th»' 
algorithms  previously  described  iii  Sections  2  and  3- 

0.  Read  Initial  data. 

1  Ghoosa  two  lint'  cities.,  x  and  y,  at  random  from  among  the  list  of 
possible  cities-, 

2.  Sliminate  the  ilxik  cities  from  thr  list  .-.  Than  put  the  remaining  ei  iias 
in  random  order 

3.  Usiiig  the  link  citia.i..  x  and  y.  as  the  initial  acyclic  permutation  on 
2  cities,  use  the  algorithm  of  Sec,ti<">n  2  to  construct  a  permutation  on  n  cities 
and  compute  the  length  of  the  tout 

i».  Go  back  to  2  and  •‘.onstruct  a  new  pemutation  Compare  each  time  wiih 
the  previous  shortest  lour  found,  saving  the  best  one  observed  in  k  trial.s, 
where  the  value  of  k  i-s  ,-ead  in  as  data,  or  else  k  is  chosen  to  be  a  mul¬ 
tiple  of  n.  For  instance,  k  -  50  or  75  or  k  2n  are  typical  values  we 
used 

5  Frint  the  be'st  tour  found  in  k  loops.. 

6.  Find  a  <;St-y  on  the  d'amete-r  cf  the  best  tour  as  follows:  Choose  any 
city,  a,  at  random;  find  the  city,  b,  on  the  iour  t.hat  is  farthest  from 
at  then  find  the  city,  c,  on  the  tour  that  is  farthest  from  b..  Select  c 
as  the  diameter  city..  ^Thla  process  could  be  repeated  more  often,  but  would 
cycle,.  The  net  effect  la  to  get  a  city  on  the  tour  that  is  '‘far  away'*  from 
the  "center"  of  the  problem. 1 

7.  Define  a  convex  subprcblem.  To  describe  how  we  do  this,  let  us  assume 

^  Vi 

that  diameter  city  c  ds  also  i^^^  the  10  city  on  the  best  tour  found  in  4 


now  CO the  ciSdtan:^  d  ironi  Lo  1^.  K^xt  we  compute  the  dlsteri'  • 

d  ftw-m  ,  i.u  t  aad  see  whether  It.  is  greatfi*  than  or  equal  to  d-^  In 
-  1-^. 

•in>  .  v.f  rect;<;e  <1  by  d'  arid  go  ono  ^lext  ws  compute  the  distance  d" 
from  j.|  ,  u:  .1.^^  ani  laako  the  asmr  wompariaor.  replacing  d  by  this  diatann' , 

■»'u  i-on<  Ir  u;  th3.s  prou'^ss”";  aa  >ong  as  d  <>ontj.n<se8  to  In'sreaseo  Onne  d  st .'ii' 
to  dccreis*  .  we  continue  the  pi-ocess  as  long  as  d  continues  to  decrease,  atii 
termi:^atfr  it  as  soon  as  d  starts  to  increase  again-.  Thus,  if  the  optima) 
tour  coneiats  of  points  on  a  '.ir^-Xa  we  would  include  tha  whole  problem  as  a 
subprobiem  Bu'  if  '  he  optimal  tour  were  in  the  shape  of  an  hour  glaasj,  as 
in  Flguin  i, .  k«  woul'i.  cut  of!  ouf*  i.obe  of  i  ^  t.tc-  13  Lustrations  of  subprofc- 
defjntd  by  ths  program  are  shown  3n  libchiblts  C,  D,  and  K.  The  raadar  wD) 
note  In  these  figure;}  that  ttie  rpsul.t:aig  subproblems  ar'e  not  always  convex  In 
the  strict  inatnr-niatica.l  senna  However,  the  resvilting  subproblems  are  auf~ 
fi-isntii  i'imp.  0  for  our  algorithm  to  solve  aasily,-  whi-^h  is  a).!  that  J.b  de.3ire.d 
c? ,  Oia  e  th^  tonvsx  subproblem  is  defin.- i»  we  print  out  the?  partial  answer  . 
If  the  ti:n''‘sx  subprcblam  includes  all  the  of  the  original  problem  we 

stop  the  pro?.«aurt  If  the  subproblem  does  not  include  all  the  oities  on  the 
liut,  W'.-  bi't  up  a  new  problem  that  consists  of  all  tha  udties  not  included  In 
tm  i-on''.fX  prodjem  pranted  out;  de-termine  t.h*<  JinV  ftit-ias  x  and  y  (being 
•'areiuJ  cfcout  r':u-  orhi-r  of  these  '-ities  in  the  nonsyrametric  case)  and  go  bavic 
to  Sjnes  a  f’initn  number  of  citisa  are  removad  each  time  tha  cities  of  a 
3ubprob.i-t'  ate  .deleted  from  the  list  of  remaining  cities,  this  process  will 
stop  .ii'.'tr  a  fini.ie  number  of  steps- 

in  another  version  of  the  program  we  employed  "double-cutting,"  that  is, 
rvi'-oving  B  aub;rctl'-‘ra  from  both  ends  of  ths  diameter  of  the  partial  tour  chosen 
Jp  6.  'i'i  _i,-i  procers  worked  reasonably  well  but  necessitated  a  longer  learoing 
loop  in  i~  to  be  certain  that  both  ends  of  the  periphery  of  the  tour  were  correct 


•  - 

There  waa  a  net  savjxig  of  rompatationa?  lljne  ,  rtowever 

Ifc  should  be  obvious  to  the  reader  that  the  pro^rm  we  l-iav*5  out  lined  was 
devised  by  looking  at  geopietrical#  henc;  s^Tomftric  and  eu-olidean,  p^'«biear. , 
’ieverthelass,  the  proixan  is  entirely  arithnjc-fci  raJ  in  nature,  and  the  Imited 
experience  we  have  rad  with  perhaps  S  a'rrisyciiwtri.'-;  problem#  andirated  ih’it  it 
works  for  these  probi'emn  ns  weli„  K-jw-j-nr.  'as  latte.!-  '-c.n'/luslon  le  highly 
tentative  and  netvds  a  gono  dsal  of  lurlher  u-.udy  and  xper liar ntat ion  bofore  it 
ean  be  firmly  shnerted.. 

uor  exparience  with  the  varloas  3iampie  probldcis  is  shown  in  Figure  5  •  The 
actual  cuts?  mad#  by  tha  manhinu  are  sfiowri  da'shad  in  Exhibits  A-F.  In  Figur-o  5 
we  have  indicated  tha-  experience  both  witii  asri  without  the  double -cut  ting  featurei. 

Of  these  problems  the  5  and  lu  city  problems  wjre  '■-omplataiy  trivial  to 
eoives  The  33  city  problem  was  interesting  in  that  the  program  fwith  aingla  cut¬ 
ting;,  had  mistakes  In  its  initial  guessas  that  ware  not  corractod  u;itil  the:  last 
cut  was  madou  In  the  42  city  problem  (with  single-cuitting) ,  the  program  actually 
obtained  the  correct  answor  after  the  first  cut,  so  that  doubtless  a  shorter 
learning  loop  would  also  have  produced  the  correct  answer u  liote  that  double- 
cutting  almost  halved  the  time  needed  to  solve  this  problem.  The  answers  shown 
to  these  four  problaros  have  been  proved  to  be  the  optimal  answare  by  various 
people  iZf  BK 

By  far  the  most  difficult  problem  is  the  57  city  problems  the  data  for 
which  was  obtained  from  the  Kand-Mctally  i?62  road  atia#  of  the  United  States. 

On  some  initial  runs  of  the  problem  we  obtained  a  tour  that  had  length 
miles  shown  in  iSSchibit  The  best  answer  that  we  found  is  that  shown  in  .:oc- 
hibit  E  whith  has  length  12, >85  miles.  Although  these  two  tours  agree  p-’«tly 
wall  around  the  periphery,  they  are  quite  different  in  the  raiddiewest.  An  an¬ 
swer  to  this  problem  that  is  mile*  shorter  has  been  obtained  by  Gordon  and  3 


•  l-Z,^ 
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Th«  answer  wad  tha  same  as  in  exhibit  E  exf^epl  for  the  mistake  in  the 
Northeast.,  dis'iusaed  in  Section  5 

Figure  ^ 


Reiter  (private  sonmiunication).  In  spite  of  repeated  attempts  our  program 
has  never  produced  their  solution^  the  probable  reasons  for  this  are  dis^- 

cussed  in  Section  6 .  It  should  be  noted  that  their  pi*ocedure  requires 
sevsra).  days  of  computation  t-o  obtain  suish  an  answer#  so  that  in  tertns  of 
eost  of  computation  the  answer  in  iixhjbit  E  is  still  probably  preferable » 
In  addition  to  the  sxperienoe  already  discussed  we  have  made  several 
runs  on  larger  problems  with  randomly  generated  data  of  dimensions  as  high 
as  90  X  90 ,  Although  the  program  works  for  such  problems  we  have  no  way 
of  comparing  how  good  the  answerii  are  with  any  other  standard  so  that  we 
shaii'  not  report  on  such  experience  her»c 
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One  method  for  >- becking  the  optSc/u;.l^y  of  proposed  answers  to  irao'elling 
salesman  prcbAaoit*  Imis  h«)?n  given  by  banr^jg.  Puikerson,  and  dohnann  Taa-^ 

method  of  -.oiurso,  ai'P-ll'/able  to  ehe  anr-vsra  whi^n  our  prograoi  givau. 

We  prop'jae  a  probabilsati .  method  tor  vhe*;kirig  on  tht#  acoara<'/  of  fhe 

result'.  I'ur  mode!  ar  »,he  i;  a>rf'"  hinomia’.  model  m  whj»di  then'  are  tw.. 

events  au-.cena  and  !i.  Fjg.x-..<  t-  W'-  tavt-,  liii'.yd  probabiiltla^  of 

obaerving  auorr  js*; 'i  in  u  tTiatt  wher't  r  .  >  t.'r*.-h4bi...ity  of  d  sun-iess  .^ach 

time  t. 
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Figure  6 


From  tha  table  it  can  be  .seen  that  events  of  probability  i  or  greater  are  ea* 
aentiaiiy  serf-aln  of  being  cbserved  in  triaiti  and  events  of  probability  .u5 
or  more  are  easentlally  tr.ertain  of  being  obaarvo-i  in  Jod  trials,  our  own  philoao- 
phy  is  to  delete  the  word  "essentially'’  and  regard  these  events  as  certain-  Th.4» 
ie  an  approximation,  and  our  aolutions  are  approximate  in  thl.B  probabilistic  sens®, 
Nevertheless,,  there  are  several  ways  of  improving  on  the  CiOnfidence  that  one  feels 
in  the  answer  so  obtained.  We  list  soma  of  these  mathods  next. 

(a’  Run  Code  .i  severaJ  times  to  get  different  ways  of  cutting  the  problem  up 
and  different  answars-.  Take  the  best  answer  of  these. 


(b  !  Gub  the  best  answer  ajariually  Jxito  sabproblaos  that  ara  dlffei’ant  from 
th«  ore's  that  Code  2  used.  Apply  Cod#  J  to  th«  aubproblems . 

fi  j  Take  adjacent  3ufcprobj.t-mwv  l.e.  -  oiiaa  that  have  a  sonuaon  ijuik  arid  wnrk 
pairs  of  thorn  using  Code  I  to  seo  if  chr  iubprobiams  define  the  same  answer  ai« 
that  giv«rn  by  Code  I. 

fd^i  From  Kigur««  hi':  pofsibie-  to  ^^tjjaato  the  probability  of  the  various 

subprobiems  bau-ig  w-jrked  -'orreotly.  By  taking  the  product  of  these  probabilities 
it  is  possible  to  estiraat-*  th«  probabiJity  of  having  found  the  tompleta  toio*  at 
random  using  Code  i,  and  thus  get  an  idea  of  how  mush  has  been  "laarned”  from  the 
puttint.  anti  factoring  procedures  of  Code  2.  The  probability  thus  obtained  is  a 
good  upper  astimata  of  the  probability  of  finding  the  optimum  tours,  and  can  be  used 
to  estimate  the  coat  of  getting  a  better  tour  than  the  one  found  so  far^ 

{e^  Of  course,  it  is  obvious  that  using  long  learning  loops  in  Code  Z  and  re-* 
peating  it  a  number  of  times  will  improve  the  reliability  estimate  that  one  van 
put  in  the  final  best  answer  found  by  the  algoritta. 

Vie  do  not  claim  that  our  program  is  infallible,  but  rather  that  it  gives  good 
answers  in  a  computationally  feasible  amount  of  computer  time.  For  iiistancssi  a 
better  tour  (yj  miles  shorter)  ia  known  for  the  57  city  problem  than  any  that  our 
programs  found.  An  explanation  of  the  failure  of  our  program  to  find  the  better 
one  may  be  found  by  examining  the  nine  sity  suoproblsm  consisting  of  the  cities 
34v  'It  6.  36,  39,  19,  3,  and  56  in  the  nor tiieas tern  part  of  the  United  States 
(see  Ibchibits  S  and  F) ,  The  optimus  tow  and  another  tour  that  is  14?  miles 
longer  are  the  two  most  probable  tours  ichosen  by  o'ur  heuristic.  The  relative 
probabilities  of  choosing  these  tours,  obtained  by  r'onning  Code  1  5C0  times 

on  this  aubproblem  are  shown  in  Figure 


-15- 


ubseirvad  ProbftbiJ,ity 


Thi*8  the.  non-optiaai  tour  i»  as  probdbla  as  the  optijnal  tour,  and  actually 

will  b«  Khosen  about  i+O  percent  of  th*  tim's:  it  happened  In  many  of  the  runs  that 
W9  made  that  the  baet  tours  found  would  contain  the  non-optiaal  tour  in  the  north- 
eastf  and  this  error  being  on  the  periphery  of  the  beet  tour  observed  caused  dif- 
ficultlesj  particularly  vihen  double-cutting  was  used.  This  tendency  of  our  pro¬ 
gram  to  choose  highly  probable  schedules*  relative  to  the  heuristic  used*  eonatl- 
tutes  a  weakness  of  the  method.  It  is  undoubtedly  for  this  reason  that  we  were 

unable  to  observe  the  Gordan4l.jitar  tour  which  was  yj  odles  shorter,  i.a.,  the 
tours  In  iSxhlbits  iS  and  F  were  probably  chosen  much  more  often  by  the  algorlthn, 
and  the  shorter  tour  is  extremely  unlikely  to  be  observed  in  the  relatival/  short 
amounts  of  computer  tine  we  used,  as  compared  to  the  times  employed  by  Gk>rdan  and 
Reiter. 

6,  PRbOF  THAT  THii  OPTIMAL  TuUK  CAh  Bi  CiioSiS;  3.  THii  SUCLIOaAi'i  GASii 

In  the  Suolldean  ease  it  ie  well-known  that  an  optimal  tour  can  never  cross 
itself.  This  follows  f^om  the  triangle  inequality,,  for  if  a  tour  is  used  that 
crosses  itself  then  it  Is  easy  to  see  how  to  pull  It  sway  at  the  crossing  point 
and  shorten  the  toiir. 

Now  consider  an  optlnal  tour  in  an  n-clty  travelling  salesman  problem 

in  the  euclidean  plane.  Vie  first  show  that  the  interior  of  the  tour  T  can  be 

n 

triangulated  by  means  of  line  se^ents  lying  completely  in  the  interior  of 
For  instance.  In  Figure  B  we  illustoate  an  eight  city  optimal  tour  that  is  so 
triangulated- 


■  lt> 


Figure  8 

The  line  sagmsints  in  tha  interior  of  tha  fcour  ara  shown  da«hedu  Note  that  the 
inside  of  the  tour  has  been  d'scfejspoaad  tntu- £i-i tingles o 

To  prove  that  this  can  always  be  done^  observe  that  it  is  vacuously  true  for 
T.  a  .and  the  only  two  possible  ^asas  for  T,  are  shown  in  Figure  9  Now  assume 
that  all  tours  ^  on  H'I  cltias  can  ba  so  triangulatad.  Consider  an  optimal 


tour  T  for  n  cities ,  Choose  a  vertex,  call  it  i, ,  at  which  the  angle 
n  1 

fomed  by  the  two  edges  of  the  tour  that  meet  at  and  which  points  into  tha 

Interior  of  tha  tour,  is  less  than  160  „  Consider  tha  cltias  and  ^ 

which  are  adjacent  to  i,  in  the  optimal  tour,  as  shown  in  Figure  10  •  We  have 

two  cases:  (a)  the  line  segoant  Joining  1^  and  ^  lies  entirely  inside  the 

tour}  or  (b)  tha  line  segnent  does  not  lie  entirely  within  the  tour.  In  ease 

(a)  we  include  the  line  sagiaent  from  i.  to  1  , 

2  n-1 


as  part  of  the  triangulation, 


-1?^ 


f-h&i/- 

and  now  the  tour^bypassss  by  making  us'a  of  this  jJLne  aegpient  consists  of 


n>l  oltlsSo  This  tour  must  be  optimal  for  that  problem  or  else  the  original 
tour  on  n  sitiee  could  be  shortened.  Hence  by  the  Induction  assumption  the 
rsmaining  tour  can  be  triangulated «  giving  a  triangulation  of  In  case 

(b)y  there  must  be  a  city,  call  it  1^,  in  the  interior  of  the  triangle  formed 
by  and  We  now  consider  the  line  seffnent  from  to  1^. 

Sither  it  lies  completely  inside  or  else  there  is  another  city  i^  in  side 
the  triemgle  formed  by  1^,  Continuing  in  this  way  we  eventually 

find  (since  there  are  only  a  finite  number  of  cities)  a  line  segment  from  i^ 
to  some  other  city  on  the  tour,  say  and  which  lies  entirely  inside  the  opti¬ 
mal  tour  T^.  This  line  segment  divides  into  two  sub-' tours  each  iiivolving 
fewer  than  n  cities.  Hence  by  the  induction  assumption  these  two  sub-tours 
can  be  triangulated,  arwl  these  give  a  triangulation  of 

By  the  same  kind  of  inductive  argument,  it  can  be  shown  that  n-3  interior 
se^uents  will  be  needed  to  perform  the  triangulation.  It  can  also  be  shown  tt^iat 
there  Is  at  least  one  vertex  that  has  no  interior  segaant  connected  to  it,  such 


as  vertex  1  in  Figure  8o 

To  demonstrate  that  there  is  a  random  order  of  the  cities  that  will  make 
the  heuristic  produce  the  optimal  tour„  we  consider  a  triangulation  of  the  opti¬ 
mal  tour  T^.  Eemamber  the  cities  so  that  city  1  is  a  vertex  having  no  interior 
segment  of  the  triangulation  connected  to  It  fas  in  Figure  8),  and  assume  that 
the  other  cities  are  numbered  in  order  around  the  to?ir  Then,  necessarily  be¬ 
cause  we  have  a  triangulatioua  cities  Z  and  n  will  be  connected  by  an,  interior 
segment  ik  and  6  in  Figure  8.)..  Hence  we  make  cities  1  and  8  be  the  first  two 
cities  considered  by  the  algorithm,  and  city  2  the  next  one  on  the  list.  The 
sequent  from  2  to  n  will  be  the  base  of  one  of  the  triangles  of  the  triaixgu- 
lationt  let  the  next  city  on  the  list  be  the  city  at  the  peak  of  that  triangle 
etc.  For  instance,  in  Figure  8  we  let  cities  1  and  8  be  the  initial  cities,  and 
introducs  the  other  cities  in  the  order  2,  3>  7*  4,  5*  and  6.  Ths  algorithm 
given  in  Section  2  will  then  produce  the  optimal  tour  1,  2,  3*  4,  b,  7>  8,  as 
the  reader  can  easily  check.  Thus  for  every  different  triangulation  of  the  opti¬ 
mal  toiir  we  get  a  different  Initial  list  that  will  produce  the  optimal  answer. 
There  are  other  initial  lists  which  will  also  produce  the  optimal  tour,  for  in¬ 
stance  the  list  1«  8,  5a  4*  4,  3)  and  7  in  Figure  6. 

At  present,  we  do  not  have  a  proof  that  for  non-euclldean  problems  there  is 
an  initial  list  which  will  produce  the  optimal  tour,  tsvertheleee,  our  experi¬ 
ence  leads  us  to  conjecture  that  such  is  possible. 
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Exhibit,  A 


5  Cmf  OPTIMAL  TOUR 
Distance  lUd 


Exhibit  B 

10  CITI  OPTIMAL  TOUR 
Distance  378 

7  8 


Exhibit  G 
List 

33  City  Problem 


1. 

Chicago,  Ill. 

17. 

Wichita,  Kan. 

2. 

Indiaxiapolis,  Ind. 

18. 

Amarillo,  Tex. 

3- 

Marion,  Ohio 

19. 

Truth  or  Consequences,  N.  Mex 

k. 

Erie ,  Penna . 

20. 

Manuelito,  N.  Mex. 

5- 

Carlisle,  Penna. 

21. 

Colorado  Springs,  Colo. 

6. 

Wana,  West  Virginia 

22. 

Butte ,  Mont . 

7- 

Wilkesboro,  N.  C. 

23. 

Lewiston,  Ida. 

8. 

Chattanooga,  Tenn. 

2^4-. 

Boise,  Idaho 

9. 

Barnwell,  S.  Car. 

25. 

Twin  Falls,  Ida. 

10. 

Bainb ridge,  Ga. 

26. 

Salt  Lake  City,  Utah 

11. 

Baton  Rouge ,  La . 

27. 

Mexican  Hat,  Utah 

12. 

Little  Rock,  Ark. 

28. 

Marble  Canyon,  Arlz. 

13. 

Kansas  City,  Mo. 

29. 

Reno,  Nev. 

l4. 

La  Crosse,  Wis. 

30. 

Lone  Pine,  Calif. 

15- 

Blunt,  S.  Dak. 

31. 

Gustine,  Calif. 

l6. 

Lincoln,  Neb. 

32. 

Redding,  Calif. 

33. 

Portland,  Ore. 

Exhibit  H 
List 

42  City  Problem 


1. 

Manchester,  N.  H. 

22. 

Denver,  Colo. 

2. 

Montepelier,  Vt. 

23. 

Cheyenne,  Wyo. 

3. 

Detroit,  Mich. 

24. 

Omaha,  Neb. 

4. 

Cleveland,  Ohio 

25. 

Des  Moines,  Iowa 

5- 

Charleston,  W.  Va. 

26. 

Kansas  City,  Mo. 

6. 

Louisville,  1^. 

27. 

Topeka,  Keuis. 

7. 

Indianapo li s ,  Ind . 

28. 

Oklahoma  City,  Okla 

8. 

Chicago,  Ill. 

29. 

Dallas ,  Tex . 

9- 

Milwaukee,  Wis. 

30. 

Little  Rock,  Ark. 

10. 

Minneapolis,  Minn. 

31. 

Memphis,  Tenn. 

11. 

Pierre,  S.  D. 

32. 

Jackson,  Miss. 

12. 

Blsmark,  N.  D. 

33. 

New  Orleans,  La. 

13. 

Helena,  Mont. 

34. 

Birmingham,  Ala. 

l4. 

Seattle,  Wash. 

35. 

Atlanta,  Ga. 

15- 

Portland,  Ore. 

36. 

Jacksonville,  Fla. 

l6. 

Boise,  Idaho 

37. 

Columbia,  S.  C. 

17- 

Salt  Lake  City,  Utah 

38. 

Raleigh,  N.  C. 

18. 

Carson  City,  Nev. 

39. 

Richmond,  Va. 

19- 

Los  Angeles,  Calif. 

40. 

Washington,  D.  C. 

20. 

Phoenix,  Ariz. 

4l. 

Boston,  Mass. 

21. 

Santa  Fe,  N.  Mex. 

42. 

Portland,  Me. 

Exhibit  I 
List 

57  City  Problem 


1. 

Akron,  Ohio 

29. 

2. 

Atlanta,  Ga. 

30. 

3- 

Baltimore,  Md. 

31. 

4. 

Birmingham,  Ala. 

32. 

5- 

Bismarck,  N.  Dak. 

33. 

6. 

Boston,  Mass. 

34. 

7- 

Buffalo,  N.  Y. 

35- 

8. 

Cheyenne,  Wyo. 

36. 

9- 

Chicago,  Ill. 

37. 

10. 

Cincinnati,  Ohio 

38. 

11. 

Cleveland,  Ohio 

39. 

12. 

Columbus,  Ohio 

40. 

13. 

Dallas ,  Tex . 

4l. 

Ik. 

Denver,  Colo. 

42. 

15. 

Des  Moines,  Iowa 

43. 

l6. 

Detroit,  Mich. 

44. 

17. 

Evansville,  Ind. 

45. 

18. 

Ft.  Wayne,  Ind. 

46. 

19. 

Harrisburg,  Pa. 

47. 

20. 

Helena,  Mont. 

48. 

21. 

Houston,  Tex. 

49. 

22. 

Indianapolis,  Ind. 

50. 

23. 

Jackson,  Miss. 

51. 

24. 

Jacksonville,  Fla. 

52. 

25. 

Jefferson  City,  Mo. 

53. 

26. 

Kansas  City,  Mo. 

5*+. 

27- 

Lansing,  Mich. 

55. 

28. 

Little  Rock,  Ark. 

56. 

57. 

Los  Angeles,  Calif. 
Louisville,  Ky. 
Memi)his,  Tenn. 
Milwaukee,  Wis. 

Mpls,  St.  Paul,  Minn 
Nashville,  Tenn. 

New  Orleans,  La. 

New  York,  N.  Y. 
Omaha,  Nebr. 

Peoria,  Ill. 
Philadelphia,  Pa. 
Phoenix,  Ariz. 
Pierre,  S.  Dak. 
Pittsburgh,  Pa. 
Portland,  Ore- 
Raleigh,  N.  C. 

St.  Louis,  Mo. 

Salt  Lake  City,  Utah 
San  Francisco,  Calif 
Seattle,  Wash. 
Spokane,  Wash. 
Springfield,  Ill. 
Springfield,  Mo. 
Tampa,  Fla. 

Toledo,  Ohio 
Topeka,  Ksuis. 

Tulsa,  Okla. 
Washington,  D.  C. 
Wichita,  Kans. 


Exhibit  J 


Data  for  the  Five  City  Problem 
Source;  Hypothetical  Example 


1 

0 

2 

30 

0 

3 

26 

24 

0 

4 

50 

40 

24  0 

5 

40 

50 

26  30 

0 

Exhibit 

K 

Data 

for 

the  Ten  City  Problem 

Source:  L.  L.  Baraehet. 
"Graphic  Solution  of  the 
Travelling  Salesman  Problem," 

0-  R-  5(1957)  841-5 

1 

0 

2 

28 

0 

3 

57 

28 

0 

4 

72 

^5 

20 

0 

5 

81 

30 

10 

0 

6 

85 

57 

28 

20 

22 

0 

7 

80 

63 

57 

72 

8l 

63  0 

8 

113 

85 

57 

^5 

4l 

28  80 

0 

9 

89 

63 

40 

20 

10 

28  89 

40 

10 

80 

63 

57 

^^5 

4l 

63  113 

80 
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EXHIBIT  N 

Data  for  the  57  City  Problee 


1  0 

t  0  Sourc:  Und-McNtllv  told  Atl«i.  38th  Edition,  Rond-McHdlly  Coepony:  1962 


3 

331 

675 

0 

4 

702 

152 

791 

0 

5 

1210 

1523 

1531 

1475 

0 

6 

655 

1074 

399 

1190 

1818 

0 

7 

210 

881 

364 

900 

1377 

447 

0 

S 

1334 

1457 

1655 

1365 

580 

1942 

1497 

0 

9 

364 

702 

685 

653 

853 

994 

527 

968 

0 

10 

227 

454 

498 

475 

1187 

‘  869 

425 

1188 

324 

0 

' 

11 

33 

695 

351 

714 

1167 

631 

186 

1311 

341 

239 

0 

12 

121 

554 

395 

584 

1205 

763 

327 

1246 

359 

106 

140 

b 

13 

1160 

620 

1422 

665 

1172 

1821 

1378 

867 

937 

953 

1254 

1054 

0 

14 

1360 

1411 

1631 

1302 

681 

1997 

1507 

101 

1013 

1172 

1366 

1236 

784 

0 

IS 

706 

899 

1027 

823 

672 

1314 

869 

628 

339 

576 

663 

701 

701 

674 

-  0 

M- 

*490. 

.7.09 

511 

730 

1139 

699 

252 

1263 

293 

255 

167 

185 

1203 

1318 

635 

0 

17 

452 

414 

723 

365 

1114 

1107 

649 

1076 

292 

225 

463 

331 

749 

1030 

505 

442 

0 

202 

60S 

556 

591 

1021 

843 

398 

1099 

175 

151 

212 

154 

997 

1149 

471 

158 

284 

0 

19 

301 

729 

80 

816 

1501 

396 

284 

1625 

655 

461 

321 

382 

1436 

1619 

997 

481 

713 

526 

0 

20 

1806 

2081 

2127 

2021 

631 

2414 

1969 

712 

1449 

1763 

1783 

1801 

1577 

811 

1240 

1735 

1720 

1617 

2097 

0 

21 

1283 

842 

1457 

666 

1414 

1856 

1481 

1109 

1091 

1056 

1295 

1162 

242 

1026 

943 

1311 

843 

1128 

1482 

18.19 

0 

22 

291 

510 

566 

487 

1032 

936 

481 

1079 

184 

109 

295 

173 

881 

1063 

467 

274 

160 

116 

556 

1629 

1011 

0 

23 

913 

402 

1036 

247 

1453 

1437 

1102 

1278 

753 

677 

916 

783 

411 

1195 

818 

924 

495 

763 

1063 

1972 

429 

647 

0 

24 

929 

316 

793 

430 

1639 

1191 

1083 

1773 

1027 

776 

960 

859 

1007 

1727 

1224 

1031 

739 

927 

865 

2429 

907 

835 

596 

0 

2S 

656 

660 

934 

586 

928 

1338 

696 

801 

393 

468 

710 

538 

562 

755 

271 

634 

298 

476 

921 

1457 

785 

365 

563 

996 

0 

26 

776 

60S 

1053 

713 

769 

1454 

1009 

652 

504 

596 

823 

658 

495 

606 

206 

775 

424 

587 

1040 

1306 

737 

485 

664 

1121 

149 

0 

27 

237 

731 

558 

721 

1069 

843 

405 

1209 

239 

270 

219 

232 

1146 

1261 

578 

84 

414 

130 

520 

1681 

1252 

246 

893 

1046 

573 

684 

0 

26 

847 

521 

1069 

393 

1193 

1488 

1045 

1049 

650 

620 

859 

735 

333 

965 

576 

871 

416 

681 

1075 

1712 

441 

565 

260 

822 

366 

404 

811 

0 

29 

2371 

2215 

2649 

2063 

1659 

3042 

2597 

1169 

2092 

2182 

2411 

2246 

1405 

1157 

1794 

2363 

2012 

2175 

2628 

1234 

1547 

2073 

1816 

2416 

1735 

1508 

2272 

1695 

30 

336 

396 

609 

37  > 

1146 

980 

536 

U70 

296 

111 

351 

218 

851 

1124 

572 

367 

125 

218 

572 

1743 

945 

114 

575 

721 

392 

518 

348 

518 

31 

717 

382 

950 

254 

1248 

1350 

906 

1111 

538 

481 

720 

596 

472 

1048 

605 

719 

277 

561 

936 

1767 

566 

445 

213 

679 

350 

459 

691 

139 

32 

452 

793 

773 

739 

771 

1060 

615 

1003 

91 

429 

429 

447 

lOlO 

1002 

356 

381 

378 

263 

743 

1367 

1164 

271 

830 

1109 

465 

550 

327 

723 

33 

769 

1060 

1090 

1037 

431 

1377 

932 

600 

412 

746 

746 

764 

953 

846 

253 

698 

687 

580 

1060 

1027 

1195 

592 

1051 

1396 

509 

458 

644 

829 

34 

513 

257 

728 

205 

1265 

1162 

718 

1206 

448 

293 

532 

398 

694 

1160 

629 

541 

160 

400 

738 

1862 

788 

297 

403 

573 

423 

554 

530 

361 

35 

1043 

518 

1145 

354 

1623 

1544 

1241 

1369 

930 

816 

1055 

928 

503 

1287 

1006 

1071 

669 

943 

1170 

2076 

358 

827 

187 

565 

760 

034 

1060 

430 

36 

487 

662 

187 

978 

1691 

216 

445 

1611 

841 

653 

507 

551 

1605 

1866 

1183 

667 

802 

712 

106 

2273 

1644 

724 

1225 

980 

1089 

1209 

714 

1273 

37 

841 

1014 

1162 

917 

584 

1449 

1004 

493 

474 

696 

818 

755 

663 

539 

135 

770 

621 

606 

1132 

1072 

905 

589 

860 

1325 

353 

204 

713 

608 

36 

451 

656 

766 

617 

858 

1102 

657 

669 

155 

323 

471 

380 

805 

910 

253 

423 

260 

242 

785 

1456 

973 

214 

661 

972 

243 

352 

369 

530 

39 

406 

772 

97 

888 

1610 

304 

366 

1730 

760 

572 

426 

470 

1524 

1785 

1102 

586 

801 

631 

105 

2202 

1554 

M3  1133 

890 

1009 

1128 

633 

1193 

40 

2011 

1833 

2289 

1678 

1500 

2682 

2237 

920 

1732 

1824 

2051 

1886 

1029 

819 

1441 

2003 

1652 

1815 

2268 

1103 

1156 

1713 

1440 

2036 

1375 

1228 

1912 

1336 

41 

1127 

1361 

1448 

1307 

207 

1735 

1290 

444 

771 

1104 

1104 

1122 

965 

522 

492 

1056 

997 

930 

1418 

739 

1207 

951 

1258 

1697 

743 

594 

1002 

990 

42 

105 

714 

230 

751 

1309 

598 

216 

1429 

459 

284 

125 

182 

1223 

1484 

801 

285 

509 

330 

200 

1901 

1326 

356 

956 

897 

720 

040 

332 

899 

43 

2492 

2662 

2813 

2572 

1317 

3100  2633 

1207 

2131 

2395 

2469 

2455 

2057 

1285 

1835 

2421 

2283 

2303 

2703 

686 

2235 

2286 

2468 

2980  2008 

1859 

2367 

2250 

44 

594 

402 

301 

559 

1794 

701 

616 

1696 

948 

521 

614 

493 

1224 

1692 

1090 

774 

671 

638 

376 

2390 

1225 

629 

006 

492 

960 

1086 

821 

897 

45 

529 

554 

606 

490 

967 

1239 

797 

906 

294 

340 

611 

411 
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55 
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1641 

1199 
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744 
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237 

303 
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56 
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38 
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492 
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1384 
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511 
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556 
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57 
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782 

1235 
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1211 
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1  2  3  *  S  6  7  8  9  10  11  12  13  16  IS  U  17  18  1«  20  21  22  23  26  2S  28  27  28 


3Sth  Idlcion,  land-NcNally  Caivaiiy;  1962 


0 

832  0 

2m  1363  0 

2916  2123  762  0 

3000  2172  879  863  0 

2712  1886  721  922  288  0 

821  101  1360  2122  2117  1829  0 

1003  221  1231  1961  2093  1803  322  0 

682  1021  2360  2831  3163  2877  1070  1069  0 

621  306  1667  2629  2291  2003  603  727  1120  0 

1169  317  1069  1811  1872  1386  376  233  1307  736  0 

1187  602  1171  1776  2062  1762  303  181  1216  908  228  0 

263  806  2113  2873  2737  2669  831  1027  963  636  1116  1222  0 

1231  669  1003  1721  1866  1371  306  268  1317  888  139  183  1268 


66  63  66  67  68  69  SO  51  92  33  36  SS  98 


